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It is proven that there exist some universal relations between the energy gap and the differences 
of the thermodynamic potential, entropy, specific heat and critical magnetic field for many two- 
and three-dimensional models of superconductivity with spin-singlet Cooper pairs forming s or d 
or g etc. states. The obtained formulae make it possible to derive thermodynamic functions and, 
in particular, the superconducting specific heat and the critical magnetic induction in the whole 
temperature range < T < T c employing the form of A(T) only. The inverse formula allowing 
us to find A(T), when the temperature dependence of the specific heat difference is known, is also 
presented. The results are referred to some obtained within the McMillan formalism, and some 
remarks on an application of the present formulae to the t-J model are given. 
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I. INTRODUCTION 

Nowadays, there exist several approaches which are able to explain properties of high-T c and other new-generation 
superconductors. In most of these approaches only one common microscopic mechanism responsible for supercon- 
ductivity in all high-T c superconductors is presented. However, this mechanism has not been univocally identified so 
far, although particularly preferred is the coupling through electron-electron interaction pj, or spin exchange [UHf]. 
as well as a modified phonon-mediated established by Anderson 0, 0, 0] , based on the concept of the resonating 
valence bond states. On the other hand, in all these approaches the pairing interaction is modeled by means of the 
exchange-mediated coupling which leads to the formation of Cooper pairs. Then,in order to explain the enhancement 
of the transition temperature, one can employ the Van Hove scenario 0, H, 0, Efl HH HE HE Hj> lla. Il6j. Recently 
it has been shown 17] that starting in the real space by setting a Hamiltonian, which contains atomic or itinerant 
electrons energy, all hopping and exchange integrals and all two-site interactions, and then solving the model Hamil- 
tonian within the mean-field approximation, one can transfer the problem into the reciprocal space. Performing a 
specific conformal transformation of the reciprocal space, one can reduce the problem to the BCS-type one, where all 
elements of the symmetry and the dispersion relation are gathered in the Jacobian of the conformal transformation. 
Such Jacobian describes a specific scalar field of the density of states which is imposed on the BCS-type isotropic 
reciprocal space, and which can contain some singularities and nodes. The Jacobian appears in all equations and its 
mean value can be identified with the density of states. Therefore, the proposed formalism is an extension of the 
Van Hove scenario and can be applied for non-s-wave symmetry systems. 

In the present paper we demonstrate that for a wide class of two- and three-dimensional models of superconductors 
with electrons or holes as carriers, and spin-singlet Cooper pairs forming a s or d or g etc. state, there exist some 
universal relations between the energy gap and the thermodynamic potential, entropy and specific heat differences 
defined between the superconducting and the normal phase. Moreover, we show that the obtained formulae allow us 
to derive the critical magnetic induction and its derivatives. We also discuss the case when strong-coupling electron- 
phonon effects are included, and we give remarks on an application of the present formulae to a relevant model of 
superconductivity i.e. the t-J model. 

II. GAP EQUATION 

Although the gap equation can be obtained in a self-consistent manner (e.g. in the Green function formalism) we 
emphasize that it defines the energy gap which minimizes the free energy, so it can be found from the free energy 
variation. Because in all known approaches the one-particle energy of eigenstates is modified by the term = 
V 7 ( e k ~~ m) 2 + I Ak| 2 , we postulate that the gap equation, which can be obtained in various approaches, formulated for 
two- or three-dimensional systems can be always written in the following general form 
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where Vk,k' is a pairing interaction, /i denotes the chemical potential and / is a certain analytic function properties 
of which are defined below. This equation should become complete with another self-consistent equation 
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which gives the chemical potential /i as a function of the filling of the conduction band n, for the normal metallic 
phase. Eq. Q after performing the conformal transformation of the momentum space [T^ | turns into the form 

4^-1 \y_ ttfi V(e-M*) 2 + |A(e^',^)i 2 v T T y '/ 

where d = 2, or 3 is the dimension of the reciprocal space, /x* defines the chemical potential measured from the Fermi 
level, so /i = ep + /i* and £ measures the energy from the Fermi level. u u and a->d mark the limits of integration 
which are the Debye energy ujjj for sufficiently wide conduction bands and one of them or both must be replaced by 
conduction band borders for partially- filled bands of width 2cj& = uj u + ujd, when ujf, < ujd- Then u> u and lu^ can be 
identified with non-symmetric cut-off parameters of the pairing interaction, and for sufficiently narrow conduction 
bands all particles get paired. When the system is almost but less than half-filling and n = 1 — 5, where S is small, 
we have 

Wd=0- — S)ub and lu u = (I + 5)uj b . (4) 

So, for some superconducting systems the Debye energy ujd must be replaced by the half band width ui, after 
including additional factors connected with the filling of the band n. Therefore, we can state that the range of the 
pairing interaction is given by the energy u> which represents ujd or u>i, according to Eq. Q and hereafter we apply 
the following notation 

LOd = nuj and oj u = (2 — n)ui 1 (5) 

where for the case of a symmetric pairing interaction as e.g. for the BCS-type models one has to put n — 1 and 
u) = u)£>. Moreover, < . . . > denotes the average over spherical angles i?' and ip' for the three-dimensional case and 
over the angle ip' only for the two-dimensional case, where §' = i 1 ) = | are consistently omitted. J(£, (f) is the 
Jacobian of the employed conformal transformation, and the pairing interaction for pure s,d 7 . . . pairing (I — 0, 2, . . .) 
can be taken in the form 

i 

v(z,?A4',<p,</) = vi(z,?) Yl yu^ra^V): (6) 

m=— I 

where Yj m (i?, <p) are the spherical harmonics. Supposing that Vi(£, £') contains only the attractive part of interaction, 
we can replace it by giVi(£)vi(£'). Then, for systems of strongly screened electrons we can assume = 1, whereas 
for a strong electronic repulsion vi (£) must be an odd function of quasiparticle energy |18l Il9l |20| . g\ is a coupling 
constant and its effective value should be obtained in a comprehensive renormalizing procedure including e.g. strong 
coupling effects and Coulomb repulsion of electrons forming a Cooper pair |2lJ . In order to be in accordance with the 
specified conditions we have to allow for the energy gap to be of the following form 

ip) = A(T) Di(-d, ip) vi(£), (7) 

where A(T) is the amplitude of the energy gap depending on temperature only, and 

i 

D0,ip)= d m (l)Y lm (#,<p) (8) 

m— — I 

expresses the spatial structure of the energy gap on equi-energy surfaces. d m (l) are normalization constants chosen 
such, that the condition < \Di{§, ip)\ 2 >= 1 is satisfied. Moreover, since for the most models under consideration the 
function / is of the form 
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as a consequence of the relation tanh(§) = 1 — 2n(x), where n(x) is the mean occupation number, we assume that in 
the limit T = it reduces to 1, whereas for A(£, i?, ip) = we have f(x, 0) ~ x when x — > 0. After employing Eq. (J7J 
and including Eq. |(HJ the gap equation J2J) can be reduced to the form 
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The obtained formula allows us to derive A(T) for a given structure of the energy gap when the Jacobian, i>z(£) and 
the function / are known. 

Let us focus now on the cases when i>z(£) = 1 for all £ or i>/(£) = const for |£ — /i* less than a few T c , then we 
denote it by vi, and when the conduction band is wide enough to imply A(0) <C Wb, and when the chemical potential 
/i is independent of temperature and equal to its value for a normal metal at T = 0, and (j,* can be omitted [22]. Then 
the above conditions allow us to ignore changes of for other values of |£ — fi*\ and leave Eq. @ out of account 
|23| . Hence, in the limit T = we have / = 1 and integrating by parts over £ the rhs of Eq. IjlUfl we obtain 
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If the last integral can be considered as constant to the postulated accuracy, Eq. Qllfl can be reduced to the form 



- = N(n) In * 
5; A(0) 

where A is positive and may be considered as independent of A(0), and 
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can be understood as a specifically defined mean value of the density of states ^3 • Because the obtained relations 
hold for an arbitrary value A(T), we have 
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Hence, subtracting Eq. (|14|l from Eq. (|10|) we can transform the gap equation to the form 
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where A = giN(n) is the dimensionless coupling parameter, F is a functional of A(T) obtained after integration over 
£ and the angles and i/3 for a given form of the Jacobian, the energy gap structure and the function /. 

Let us investigate the form of the gap equation (|10f> in the limit T = T C when A(T C ) = 0. Integrating by parts over 
£ and taking into account the above assumptions we can transform it to the following form 
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which can be written as 
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where B and 6 are positive and, due to the postulated accuracy, may be considered as independent of T c , and 

<[J((2 - n)u, , + J(-nw, 0, ^)| 2 ) 



6 = 



/ ((2 - n)f ) J((2 - n)w, ^ 0, ¥>)] |A(0, ^)l 2 ) ' 

is independent of uj. Note that the characteristic ratio 



(18) 



A(0) A (b-l\ 

exp — — , (19) 



T c B > \ A 
unlike in the BCS theory, is not a universal constant. 

III. ENERGY GAP AS A FUNCTION OF THE COUPLING PARAMETER 

The difference between the thermodynamic potential of the superconducting and the normal phase, which is equal 
to the free energy difference when the chemical potential of both phases is identical, can be derived from the relation 



AQ(T) = -N(n) [ 
Jo 



^JdA', (20) 



where we have taken into account that gi = X/N(n). The formulas l|12|) and l|17f) allow us to express A(0) and T c by 
means of an arbitrary A' varying from to A in the following forms 

A v (0) = Aexp^-l^ and T x ,=Bexp(~Y (21) 

where, to the assumed accuracy, we can also put A = A\lo and B = B\u>, where A\ and B\ are real numbers. In order 
to keep the standard notation, hereafter we assume that A>(0) = A(0) and T\ = T c , and moreover A\(T) = A(T). 
Note that T c is a decreasing function of b, and for the BCS model 6=1. Moreover, after eliminating A from these 
formulas we can find the characteristic ratio (|19f) as a function of lu. Then 



b—l 



6-1 



^)=W(^r or M =AS W« Y , ( 22 ) 
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and it becomes an increasing or a decreasing function of lo when 6 > 1 or 6 < 1, respectively. 

Employing the formulas (|21|) we can find the thermodynamic potential difference in the limit T = 0, and we state 
that it can be expressed in the following universal form 

Aft(0) = -^N(n)A 2 (0). (23) 

In order to derive the thermodynamic potential difference for non-zero temperature we have to express A\r(T), i.e. 
the solution of the gap equation (|15fl for A' varying from to A, by means of A(T') where T" can change from to 
T c . Let us suppose that we can find the energy gap for all values A' as a function of reduced temperature T/A^'(0). 
Then introducing the reduced gap in the form 

T 



A\i 

Ay(0)J __ A v (0) ' 
where, moreover, A\(. . .) = A(. . .), and employing Eq. (|21(l we can transform Eq. Ijl5|l to the form 
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Since this equation has the same form for all values A' from to A when T, u> and n are fixed, all forms of A\> (T / Ay (0)) 
derived from Eq. (|25[) must be equivalent. They become identical for temperatures normalized by corresponding critical 
temperatures denoted as r. This observation allows us to express the reduced gap derived for an arbitrary A' by means 
of the reduced gap A, which represents the real energy gap of a superconductor 



Aa'(G) 



(26) 



Hence, after some transformations we find that the value of an arbitrary reduced gap Av at T/Ty can be expressed 
by means of the function A as 



A v 



/ T 



A 
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To 
Ty 



Av(0) 
A(0) 



and the inclusion of Eqs. I|21|) and (|24[1 results in the following fundamental relation 

exp 



^' (|:)^ A (|- exp 



(27) 



(28) 



Thus, in order to find the values of the energy gap for A' < A and for a fixed temperature T, we have to know the 
real energy gap values for all higher temperature values if b > i, and for all lower temperature values if b < \. For 
b = \ the temperature dependence coincides. 



IV. UNIVERSAL RELATIONS 



The derivation of the relation (|28|l creates new possibilities to find the thermodynamic potential difference and its 
derivatives with respect to temperature as functionals of A(T). Substituting Eq. I|28|l into formula l|2U|l . after some 
calculations, we find the thermodynamic potential for b ^ ^ in the form 



Afi(T) = - 



1 



26-1 



N(n)T^ 



A 2 (T'^ 



-dT\ 



t (V) 



where T* = T c if b > \ because A(T) = for T > T c , and T* = if b < §, whereas for b = \ we have 

An(T) = -±N(n)A 2 (T), 



(29) 



(30) 



and hereafter we do not consider this case. The formula (|29|l can be also written in an another quite equivalent form 



Afi(T) = 
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N(n)T~ 



A 2 (T') - A 2 (T) 



1 



A^ dT > - -N{n)A 2 {T) 



l)^e(b-l 



(31) 



which can turn out to be more convenient for investigations. Q(x) is the Heaviside step function. Since 

AS(T) = -AaO(T) and AC(T) = -T-=Ml(T), (32) 



where AS(T) and AC(T) are differences of the entropy and the specific heat, respectively, differentiating the formula 
(I29"|> we obtain 



AS(T) 



-NlnjT^r- 1 
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which can be also written in an another quite equivalent form 
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Now, employing formula ilM-Sl we can find the specific heat difference in the form 

v ; (26 - 1) 3 v ; J T ( T /) 25^+1 (26 — l) 2 T 2b-l y ' K 'dT y h K ' 

which can be also written in another quite equivalent form 

= 2(3-26) ^ / T ' A 2 (r)-A 2 (D 3-26 A^T) / ^ / _ 1 

- - iV(rc) A (T) A (T) . (36) 



26-1 v ' x ' dT 

The obtained formulas include the parameter 6 which is an arbitrary positive number. Hence, if we could identify 6 
we would be able to find the forms of the thermodynamic potential (free energy), the entropy and the specific heat 
for the superconducting phase based upon A(T) shape. Note that for BCS-like models 6=1, and this case has been 
considered by us in details before [UlU 127| . On the other hand, if we knew A(T) and AC(T), then by virtue 
of the presented formalism we could estimate 6. We emphasize that in general N(n), i.e. the density of states defined 
above for the superconducting phase, is different from that one defined for the normal phase. This difference can be 
explained as a modification of the effective mass, and must be included in all normal-phase thermodynamic functions 
when we derive the thermodynamic functions for the superconducting phase from the above formulas. This procedure 
is necessary in some approaches 0,113 to eliminate non-physical results, e.g. negative values of the superconducting 
specific heat. Note that the obtained Eqs. (|33[l - i|3()[) imply a formal relation between AS(T) and AC(T), which can 
be written in the form 

AC{T) = 2T~i AS(T) 2b^i N{n) ~tk A2{n (37) 

In order to consider the specific heat in low temperature limit (T — > 0), we have to take into account the normal 
specific heat Cn(T) which for this case in the frame of the presented formalism can be found in the form 

/>oo 

C N (T)=T I N(2Tu)u 2 cosh~ 2 (u)du, (38) 

where 



o 



N(2Tu) = (J(2Tu, i?, if) + J(-2Tu, >&, ip)) (39) 

is another form of the density of states (cf. Eq. I|13() '). Since N(2Tu) depends on T, Cn(T) it can become a 
complicated function of temperature. According to Eq. (|3TJl and including Eq. (|3*8|l the superconducting specific heat 
can be expressed in the form 

Cs(T) = 2^ AS{T) 26^T iY(n)A(0) ^ A(T) + Cn{T) - (40) 

Hence, we state that the formula for the superconducting specific heat can contain several different terms which 
determine the temperature dependence. For the BCS model with a s- or d-pairing the superconducting specific heat 
is proportional to dA(T)/dT only, since for 6=1 the term AS'(T) is proportional to T and it is compensated by 

C N (T) 00123 

Let us now consider the specific heat jump at the critical temperature T = T c . In order to do that comprehensively, 
we have to assume that for T — > T c the energy gap A(T) — » A(T C ) and A(T C ) = for the second-order (continuous) 
transition or A(T C ) > for the first-order (discontinuous) transition, and for T > T c the energy gap vanishes, i.e. 
A(T) = 0. Then, the specific heat jump can be found from Eq. (|37(l in the limit T — > T c as 

AC(Tc) = 2^1 AS{Tc) ^ 2lhi N{n) T lhn T c A ^ A ^ ( 41 ) 



Since for continuous phase transitions AS(T C ) = 0, the obtained formula coincides with that known for BCS-type 
models 
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V. RECONSTRUCTION OF THE ENERGY GAP 

In the previous section we have shown that the thermodynamic potential, the entropy and the specific heat differ- 
ences can be defined based upon the temperature dependence of the energy gap only. In this section we show that also 
the energy gap A(T) can be reconstructed from the specific heat difference if the values of 6 and AS(0) are known. 
Employing Eqs. (|3~2"|) and (|3Tjl . after some transformations, we can express the energy gap by means of AC(T) and 
AS(0) in the form 



A 2 (T) = A 2 (0) + 5-^TAS(0) + -l 

N(n) N(n) J 



f 


(3 - 26) / 


Jo 


Jo 



r 



AC(T") 



dT" - (26 - l)AC(T') 



dT'. 



(42) 



Hence, for the cases when Cjv(T) ~ T q , where according to the Nernst law q > 0, Eq. I|42|l can be reduced to the form 



A 2 (T) - A 2 (0) + i-^TA5(0) - 3 ™ f b ^ )q TC N {T) 



N{n) 



N(n)q(q + 1) 



+ 



N(n) Jo 



r 


(3 - 26) / 
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T' 



C S (T") 



T 



—LdT"-(2b-l)C s (T^ 



(43) 



which for BCS-type models, when A5(0) = and b = q = 1, includes Cs(T) only [2f 



VI. CRITICAL MAGNETIC INDUCTION 



The formulas obtained in Section llVl allow us to find the critical magnetic field H C (T) or the critical magnetic 
induction B C (T) = [1qH c (T) and their derivatives for superconducting systems as functionals of A(T), only. In 
accordance to the thermodynamic relation 



G S (T,H = 0) - G N (T,H = 0) = -i Mo iJ 2 (T) 



(44) 



after taking into account that Gs(T,H = 0) — Gn(T,H = 0) = Afi(T) when the chemical potential is independent 
of temperature and identical in both phases, the critical magnetic induction can be presented in the form 



B c {T)=B c {0)\ 



1 -2AQ(T) 
7V(n)A 2 (0)' 



(45) 



where 



B c (0) = V^oN(n) A(0), 



(46) 



and Af2(T) should be taken in one of the forms (|31|I - H33|) . In a similar way employing standard relations and including 
forms of Eqs. H35fl - H38[) we can express derivatives of the critical magnetic induction in the form 



dB c (T) 
dT 

d 2 B c (T) 
d 2 T 



AS(T) 



2 y/-AQ(Ty 
AC(T) 



(A5(T)) 2 
ly-An(T) (-AO(T))i 



(47) 
(48) 



which should be useful to estimate values of the introduced parameter 6 for some unconventional superconductors by 
means of experimental data. 



VII. REMARKS ON MCMILLAN'S CASE 



In the formalism developed by McMillan 21], where strong-coupling electron-phonon effects were included, the 
critical temperature was derived in the form 



T c ~ w D exp 



1 + A 
A-m*(1 + A) 



(49) 
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where the parameter fi* expresses the Coulomb repulsion between electrons forming a Cooper pair, and for known 
materials /i* ~ 0.13, cup is the Debye energy which is the actual cut-off parameter for conventional superconductors. 
Note that the critical temperature T c = for A = ^ > 0, so in obedience to the presented consideration we should 
rather use an effective coupling constant which varies form to its real value when the coupling interaction is included. 
Such effective coupling constant is A fr = a" 1 [A — fi*(l + A)]. Hence, putting it into Eq. 1)49(1 we obtain 

T c ^^exp(- T ^)exp(--A-), (50) 

where 6 _1 = a(l — /!*) > and a is positive constant which ensures that A(0) keeps the form 1(21(1 . Hence, the 
formalism developed by us can be applied to the McMillan- type superconductors. 



VIII. APPLICATION OF DERIVED RELATIONS TO HIGH-T C SUPERCONDUCTORS 

The mechanism of superconductivity for HTSC based on antiferromagnetic correlations was firstly proposed to 
predict antiferromagnetic behaviour for pure La2CuC>4. As soon as these materials are doped, giving a few per cent of 
holes in the Cu02 plane, the antiferomagnetism is destroyed. As the number of carriers 5 is always small, they cannot 
completely screen the Coulomb interaction on copper site. Thus the copper electrons can move but are always subject 
to strong on-site repulsion. The introduction of additional holes causes that in sufficiently low temperatures they are 
coupled by superexchange interaction, through the oxygen atoms, forming Cooper pairs |30l 1.3 ll l32t l3.'l 13 ll l35 | . The 
similar situation occurs in the compound YBa2Cu306 which is an insulating antiferromagnetic. In this so-called t-J 
model the gap equation can be obtained in the form Q by performing a BCS linearisation of the J term. Then it 
can be expressed in the form 

2 1 ^ 2 tanhQ3£ k /2) 

-j = N^ ' (51) 

k 

where 

7k = cos(fc x a) ± cos(kya) (52) 

and ± refer to s-wave and ci-wave superconductivity, respectively. Note that 7k represents a factorable pairing 
interaction which determines the form of the energy gap. So, it is compatible with the form of the pairing interaction 
assumed in the presented approach. Moreover, Cyrot et al. 22] in order to derive the critical temperature introduced 
some restrictions on the form of the density of states separating in the integral form of the gap equation ((51(1 with 
A = two domains of integration: the first one where the density of states can be approximated by a constant value, 
and the second one where the hyperbolic tangent can be replaced by 1. Declaring that their solution is valid only for 
the range of doping S > 0.05, when the order parameter has the s symmetry, they obtained the critical temperature 
in the form 

T c = l.UW(fM) exp[-l/JiV(^)] (53) 

where JN{y) corresponds to the dimensionless coupling constant A, and W(/x) is proportional to lo c (cf. Eq. I(17|) 
and Eq. J2U). Moreover, they showed that at zero temperature ^(fi) J A(0) = 1.76T C . So, the t-J model within the 
developed approach |22j reveals BCS-like properties with 6=1. 

In order to employ the found universal relations to the t-J model we have to put 6=1, since the other quantities 
or representing them expressions are not involved in these relations. Moreover, let us emphasize that in our method 
the summation over k is replaced by the integration over £ = e — /i for an in general non-symmetric pairing interaction 
by virtue of the conformal transformation, in which the Jacobian represents the scalar field of the density of states 
[l7| . It causes that our method is more precise than the one discussed above, and it allows us to predict that b can 
be different than 1. Note that the three substantial assumptions mentioned above imply the following facts: the 
constancy of the density of states corresponds to putting the Jacobian constant, the inclusion of a symmetric pairing 
interaction is equivalent to n = 1, and for u> = lod, same as the hyperbolic tangent, the function / can be replaced by 
1. As a consequence the parameter b given by Eq. 1(18(1 reduces to 1. 



IX. CONCLUSIONS 



The presented formalism proves that there exist some universal relations between the free energy, entropy, specific 
heat differences, and the energy gap amplitude for a wide class of models of superconductivity. In our previous 
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papers 0, HE H3 we have shown that applying the formulae (|35|l or (|36|) established for the case when 6=1 
towards a general s-paired BCS case, when the magnetic field and superflow are taken into account, and for a pure 
d-paired BCS model, we obtained correct forms of the superconducting specific heat, where the linear (with respect to 
temperature) terms were entirely eliminated by Cjv(T). Therefore we expect now that these relations will be verified 
experimentally, and they will turn out to be helpful in explanation of common elements of the mechanism of low- 
and high-T c superconductivity. On the other hand, we remind that we imposed some restrictions on conditions which 
must be satisfied by the superconducting system. Hence, the presented formalism can be applied to superconducting 
systems when the cut-off parameter lud exceeds A(0) or u/A(0) by far, the chemical potential is identical in both 
phases, and the pairing interaction can be taken in a form corresponding to BCS-type models, which can be always 
attained by performing a BCS linearisation of the appropriate interaction term. 
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